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Abstract 

In a delayed coincidence experiment, for example, the recent reactor neutrino oscillation experiments, a 
precise analytic determination of the delayed coincidence signal efficiency and the accidental coincidence 
background rate is important for the high accuracy measurement of the oscillation parameters and to under- 
stand systematic uncertainties associated with fluctuations in muon rate and random background rate. In this 
work, a data model is proposed to describe the full time sequence of all possible events on the live time axis. 
The acceptance of delayed coincidence signals, the rate of accidental backgrounds and other coincidence 
possibilities are calculated by assuming that all of the 'net muons' are uniformly distributed on the live time 
axis. The intrinsic relative uncertainties in the event rates are at the 10 -5 level for all combinations. The 
model and predictions are verified with a high statistics Monte Carlo study with a set of realistic parameters. 
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I. INTRODUCTION 



The delayed coincidence method is used very broadly in nuclear and high energy physics ex- 
periments. The existence of a delayed signal greatly mitigates the critical requirement on low 
random background levels. The experimental observation of neutrino oscillation is one of the im- 
portant breakthroughs in studying the Standard Model physics in thepast half century. Recently 
three neutrino experiments, Daya Bay, RENO and Double CHOOZ yJ-yD adopted this technique 
again to measure the neutrino mixing angle 9\ 3 precisely. Those anti-electron neutrinos from re- 
actors are distinguished by detecting the positron and neutron from an inverse beta decay (IBD): 
p + v e — > e + + n. The whole signal is composed of a prompt positron signal and a delayed neutron 
capture signal, where the prompt positron has a detectable energy between 1 to 10 MeV and the 
delayed neutron spends about 29 fxs in thermalization to be captured in the gadolinium (Gd) - 
loaded scintillator region where several photons of total energy 8 MeV are released. 

In the past year, the Daya Bay Experiment [1] and the RENO Experiment [2] have reported 
their first 9 13 measurements with delayed neutron capture signals on Gd. So far the accidental 
background is one of the dominant backgrounds. They both quoted the uncertainty for this back- 
ground is at the 1% level. Their calculation methods are the same and can be found in their 

A n 

public presentations and papers, e.g. [1] and fl2fl. The formula for accidental background number, 

N accidental, is repeated here for discussion: 

^accidental ^n—like ^ [1 CXp(R e + _ Hk e X T coincidence^- 

In a low random background level environment, i.e. the number of neutron-like signal (N n _uke) 
and the rate of positron-like signal {R e +_nke) are both low, or with a very small coincidence win- 
dow (^coincidence), all the single signals will be very well separated in time on average, for which 
this simple formula is valid. But when the random background rate is high, it will become more 
likely to have other cases, like three-fold coincidence events. In this special physics context, the 
visible energy ranges of the prompt positron and the delayed neutron overlap, giving this formula 
an intrinsic bias, i.e. R e +_uu e will be higher than the real R e +_i ike by including n — like signals. 

The Double CHOOZ yD experiment estimated the accidental background by shifting the co- 
incidence window by 1000 fis, which gave a larger statistical uncertainty of about 10%. There 
is also an additional systematic uncertainty in this off-window method: because of the spallation 
products produced by muons, the accidental background rate is lower when further in time after 
the previous muon, i.e. in the off-window. 



For these three experiments, the Gd-loaded liquid scintillator region is enclosed by a pure liquid 
scintillator region which has an almost equivalent mass as the Gd region. There are proposals to 
use the hydrogen neutron capture signal to further enhance the neutrino signal statistics, and to 
even aim at the comparable precision level as the Gd capture result. The effect mentioned above 
may only have marginal effect on the published 6> 13 measurement since the accidental background 
itself is small. The neutron capture on hydrogen will release a 2.2 MeV photon and the average 
capture time is around 200 fis. Apparently the delayed neutron signal energy is significantly lower 
than that in Gd-loaded scintillator which would allow more neutron-like singles. In addition, 
the coincident window chosen will be much longer. Correspondingly the accidental background 
will be much higher. This necessitates a precise analytic delayed coincident signal efficiency and 
accidental coincident background rate calculation. 

Another interesting issue is that the full data-taking time is chopped into slices by cosmic 
muons. Some complex muon veto cuts are used in these analyses, i.e. IBD events are rejected 
if they tightly follow a muon. The length of veto time relies on the amount of the muon energy 
deposited or the location of the muon hits. It is imaginable that while the muon rate is high, there 
is no acceptance for a delayed coincident signal or an accidental coincident background. But how 
exactly the muon rate is going to affect the signal selection etc. is not transparent. 

In this paper, we propose a data model to depict the full time sequence of all kinds of possible 
events in Sec. HH and with a possible signal selection scheme, the acceptance of delayed coinci- 
dence signal, the rate of accidental background and other coincidence possibilities are calculated 
analytically in Sec. Ull In Sec. [IV] a high statistics simulation is done to verify the model and 
predictions. 

II. DATA MODEL AND SIMPLIFICATION 
A. Data model 

A model is built up to demonstrate different types of events. Although the context used is closer 
to the physics of the reactor neutrino experiments, there is no loss of generality. In the following 
we use e and n to represent the positron and neutron from an IBD reaction respectively, s for 
random background (or singles) and ji for muon. The features of the event types are described 
below: 
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1) The most interesting event is the delayed coincident event including a prompt signal, e, and 
a delayed signal, n. The event rate of IBDs is assumed to be Ribd- The time between the prompt 
and the delay signal can be complex. For the three neutrino experiments mentioned above, the 
neutron capture cross section varies with its momentum, so that an exponential distribution is only 
a pool approximation for the neutron capture time. Here an abstract form Pdeiay(t) is used to 
represent the delayed signals' time distribution. 

2) Random signals, or singles, include those decays from residual radioactive nuclei in the 
detector or other non-correlated detector noise. Random signals occur with a uniform distribution, 
or equivalently the time interval between two random signals follows an exponential distribution 
with an average value of 1/R S , where R s is the singles rate. 

3) Cosmic Ray muons are always present. They themselves are usually easy to identify because 
of their high energy deposit. The spallation products of muons, however, can cause a brief increase 
of random signal rate, which causes a trouble for the delayed coincidence event searching. There- 
fore a passive veto window is applied for each muon. The chance to produce spallation products 
depends on how far a muon track is from the sensitive region and the energy deposit. Based on 
this knowledge, different muon veto window lengths can be chosen to maximize the live time and 
at the same time have a good control of spallation background. 

These different types of signals in the full data-taking time axis are shown in Fig. \T\ First we 
discuss two muon events, A and B, on the full time axis. Muon A is supposed to be closer to the 
sensitive region of the detector than muon B, and so a longer veto window is applied to muon A. 
A pair of delayed coincidence signals (a positron and a neutron) occurs between muon A and B, 
where T c is the coincidence window used to search for delayed coincidence events and is usually 
at the scale of the average arrival time of the delayed signal. To ensure the spallation background 
of muons will not contaminate the search for the prompt signal as well as the delayed signal, the 
time between a prompt signal and the next muon event cannot be smaller than T c . The total dead 
time introduced by a muon is the veto time plus T c . After muon B, a single signal occurs, since no 
delayed signal is found within T c after the single signal. Finally we show a case where two muons, 
C and D, are very close to each other in time such that their dead times overlap. In this example 
all coincident signals and singles are all well separated. We will discuss situations of overlap in 
the following sections. 
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FIG. 1 : (Color online) A graph to demonstrate the full model of three types of signals and how they can be 
simplified for an analytic calculation. See the text in Sec.|II]for detailed explanations of the notation. 

B. Simplification 

The first attempt at simplification is to project all signals from the full time axis to the live time 
axis as illustrated in Fig. \T\ A muon with its dead time is abstracted as a virtual point on the live 
time axis and is counted as one net muon, like muons A and B in the plot. Muons C and D have 
overlapping dead times, so they together are abstracted as one point on the live time axis and only 
counted as one net muon. The other signals are all transparently moved to the live time axis. Since 
every muon dead time with a buffer window T c is already removed from the live time axis, the full 
length of the live time axis is the total effective time of the prompt signal searching. 

The important assumption making this analytic calculation possible is to postulate that the net 
muons are uniformly distributed on the live time axis and the time interval follows an exponential 
distribution with an average value of where is the net muon rate on the live time axis. 

In the calculation below, we always have Ribd <S R s - For the three reactor neutrino experi- 
ments we took as examples, R s is about 4 orders of magnitude higher than Ribd- 



C. Searching for delayed coincidence events and other backgrounds 

In this discussion we will use the simplest method to search for delayed coincidence pairs: 
as depicted in Fig. \T\ for the positron and neutron, whenever a signal appears and is not within 
a muon dead time, a window will be opened for a duration of T c during which a second signal 
will be searched for. Two random signals can form an accidental coincidence background if their 
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separation is less than T c . 

Except to predict the event rate of delayed coincidence events and accidental background, other 
types of combinations are also interesting to check. The following cases will be calculated: 

1- fold coincident events: s, e, n. 

2- fold coincident events: ss, se, sn, en, es, ns. 

3- fold coincident events: sss, sse, ses, sns, ssn, sen, ens, esn, ess, nss. 

HI. THE METHOD OF CALCULATION 

A. Two Steps of Calculation 

The calculation of these event rates is divided into two steps: a) determine the probability of 
a type of signal to start a coincidence searching window; b) determine the probability that there 
is a second or third signal in the searching window for 2-fold or 3-fold coincidences or there is 
no other signal for 1-fold events. In the following sections the starting probabilities of a random 
background signal (P s - S tart), positron (P e - S tart) and neutron (P n _ start ) will be calculated first and 
then followed by all kinds of event rates. 

B. Starting probability 

On the live time axis, each signal except a muon can start a coincident searching window, as 
long as it is not in the previous coincidence searching window. Note that on the live time axis, 
Rjbd and R s are exactly the same as on the full time axis. First let us consider how a single event 
may start a searching window in different situations. In the formulas below, is the time from a 
signal event to its previous net muon event. 

Case a) As shown in the top panel of Fig.[2l when < T c , if there is no other signal between the 
single event and the muon, a searching window will be started by the single event. The probability 
of this situation is: 



P 



a 




-(Rs+Rv)Tci 



(1) 



R s + R, 
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FIG. 2: (Color online) Three situations of how a single event starts a searching window. See the text in 
Sec.Hn|for details. 

where P(0|i? s t At ) is the probability that there is no other single event in between, which is calcu- 
lated as a Poisson distribution with a mean of R s t^ L and count k = 0, and the second term Pit^) 
gives the probability of finding a muon at before the random background being considered, 
which is calculated according to an exponential distribution with a rate of R^. 

Case b) While t M ^ T c , if there is no other signal in a T c window before the single event being 
considered, a searching window will be started, just as the panel b of Fig. |2] shows. The starting 
probability of this situation is: 

POO 

P b = P(0\R s T c )- / P(t M )-dt„ 
(R s T^ k r °° 



" e fc=0 



(2) 



where the first term corresponds to no signals in T c and the second term is to ensure t M ^ T c . 

Case c) In case b) when t At > T c , if there is another scenario that the target single can start a 
searching window as depicted in the panel c of Fig. |2] There is a random signal, si, in front of 
the target single event; however it is occupied by a previous searching window started by s2. The 
probability of this situation is: 
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i>l c poo 

P c = / dtR s e- R ^-[l-P(0\R s t)]- / Pit^-dt^ 

JO JT c +t 

= r dtR s e~^ ■ [1 - { -^e- R ^\ k=0 } ■ T R,e- R ^-dt, 

t/ * </ Tq ~\~t 

= Rs c ~RuT c ^ _ e -(fl s +R M )T c j _ fl» C ~ RuTc [l - e - { ? Ra+R ^ Tc \, (3) 

where the first integral gives the probability to find the first random background signal si at time t 
in front of the target single event, the second term calculates the probability of having at least one 
random background s2 in the early t window, that is T c away from the target single event, and the 
last integral just gives the probability of finding a muon at some time larger than T c + t. 
Finally the starting probability is, 

Ps-start = Pa + Pb + P c - (4) 

There should be higher order corrections after P c , but, with the example parameters used in the 
simulation presented later, they are estimated to be five orders of magnitude smaller than these 
leading terms. 



The starting rate is 



As for IBD signals, it is 



Rs— start Rs ' Ps— start' (5) 



e— start 



Ps— start ' (6) 



where Rjbd <^ R s and there is a detection efficiency e e for positrons. Then the rate of prompt 
signals, i.e. positrons starting a time window, is similar to the rate of random backgrounds starting 
a time window: 

Re-start = RlBD ' Ps-start ' £e- U) 

The situation of a neutron starting a searching window is a bit complicated, because by design 
it is always possible that there is a prompt positron signal ahead of it. It is only possible when the 
positron is failed to be detected. However the method to calculate P n - S tart is the same as above 
and the final result is listed in the appendix. Then the event rate with neutrons as a start is: 

Rn— start — RlBD ■ Pn— start- (8) 
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C. Construct an Event 



After the starting probability or starting rate of a type of signal is known, the rate of accidental 
background, detectable IBD pairs and other cases can be calculated. 

For example, the accidental background rate, i.e. the single-single combination in one T c , is 
just the rate of one single event to start a searching window multiplied by the probability of a 
second single event appearing in the same window: 

R(ss) = R s _ start -P(l\R s T c ), (9) 

where R(ss) is used to describe the event rate of ss-type backgrounds and P(1\R S T C ) is the Pois- 
son probability of one count with a mean of R S T C . IBD events are not explicitly required to be 
vetoed from this with R IB d ^ R s - 

The IBD event rate can be obtained in a similar fashion: 

f Tc 

R(en) = Re-start ■ £n\e \ P delay (t)dt • P(0\R s T c ) (10) 

Jo 

= RlBD ■ P S —start ' &e ' &n\e ' I Pdelay 

(t)dt-P(0\R a T c ), (11) 

Jo 

where e n \ e is the efficiency of neutron detection after a positron is detected, and the random back- 
ground is explicitly required to be outside this searching window. The complete IBD event detec- 
tion efficiency can be expressed as: 

e, BD = |^ (12) 
The results of other combinations are listed in the appendix. 



IV. VERIFICATION WITH MONTE CARLO 



Since the muon veto cut is quite complicated, a high statistics Monte Carlo simulation (1E10 
events) is done to verify the assumption and the predictions. Three types of events are produced: 
IBD, random background, muon. They are generated on the full time axis according to three 
uniform distributions. Their event rates on the full time axis are chosen to be very similar to 
those of recent reactor neutrino experiments. Usually a veto time of 400/is is applied for each 
muon, but 0.05% of them can be shower muons or closer to the sensitive region, for which a 
one-second long veto is applied. The neutron capture time distribution Pdeiay(t) is replaced with a 
simple exponential distribution with rate = A. The detection efficiencies of the prompt and delayed 
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signals are also included. All the parameters used are summarized in Tab. U The measurements 
with the simulated sample are compared to the predictions and are listed in Tab. HH where the 
discrepancies are all within a 3 a range. 

One interesting discovery is that the net muon rate on the live time axis and the real muon rate 
on the full time axis are the same within statistical uncertainty. With the analytic formulas all 
errors on singles rate and muon rate can be easily propagated to the final result. 



R s 


Total muon rate 


veto 


Shower fraction 


Shower veto 


RlBD 


e e 




T c 


1/A 


50 Hz 


200 Hz 


400 fis 


0.05% 


1 s 


0.1 Hz 


1 


0.8 


400 /xs 


200 fj& 



TABLE I: Monte Carlo simulation parameters, including singles rate, muon rate, IBD rate, IBD delayed 
signal capture time constant and signal selection cuts and efficiencies. 



Rate [Hz] 


R(s) 


R(ss) 


R(se) 


R(sn) 


R(sss) 


R(ens) 
+R(esn) 


R(sse) 
+R(ses) 


R(sns) 
+R(ssn) 


Mea. 


48.0833 


0.96165 


0.0010460 


0.0001108 


0.009626 


0.001342 


2.10E-5 


2.27E-6 


Sta. Err. 


0.0015 


0.00022 


7.2E-6 


2.3E-6 


2.2E-5 


1.1E-5 


1.4E-6 


4.7E-7 


Pred. 


48.0856 


0.96171 


0.0010499 


0.0001093 


0.009617 


0.001330 


2.10E-5 


2.19E-6 


Rate [Hz] 


R(e) 


R(en) 


R(es) 


R(sen) 


R(ess) 


R(n) 


R(ns) 


R(nss) 


Mea. 


0.029657 


0.066447 


0.0005991 


0.0008715 


6.91E-6 


0.012719 


0.0002616 


3.11E-6 


Sta. Err. 


0.000038 


0.000057 


5.4E-6 


6.6E-6 


5.8E-7 


2.5E-5 


3.6E-6 


3.9E-7 


Pred. 


0.029647 


0.066525 


0.0005930 


0.0008735 


5.93E-6 


0.012770 


0.0002554 


2.55E-6 



TABLE II: Measured event rates of simulation sample and their predictions. For each type of combination, 
Mea. gives the measured result, Sta. Err. is its statistical error, and Pred. is the predicted value, sse and ses 
are not distinguished, as well as sns and ssn. 



V. CONCLUSION 

We brought up a complete mathematical model for the signals distribution on the full time 
axis for delayed coincidence experiments, and further projected this complicated model onto the 
live time axis. An analytic calculation was done by assuming all the net muons are uniformly 
distributed on the live time axis. The intrinsic relative uncertainties for all combinations' rates 
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are at the 1CT 5 level. The predictions were verified with a high statistics simulation with a set of 
realistic parameters from the Daya Bay, RENO and Double CHOOZ experiments. The method 
proposed here is precise and can be used in these reactor neutrino experiments and other related 
experiments. In the appendix the event rate formulas for all the cases are listed. 
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APPENDIX 



All 1-fold, 2-fold and 3-fold coincident event rate calculation formulas are listed in Tab. HTH 
where the first column indicates the combination concerned and the second column shows the 
expression for calculation. P(k\m) is used to describe a Poisson distribution with a mean of m 
and a count of k. The delayed signal follows a simple exponential distribution with time constant 
A. 
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^~^s — start 


R, ■{„*£„ [1 e-C^+^cl + e -(R s +R„)T 0+ 

Rs+R,, C l L C J 2R S +R„ e l L e \i 


Re— start 


Rs r -R,,T r u -(Rx+R^Trl R 3 r -R,,T r u -(2R«+R„,)T n u 
Rs+R„ C \- L C J 2R a +R [l e l L e \S 


Rn— start 


RiBD{e n R^[{l e e |nj Rs+Rfl + e e \ n Rs+R ^ +X \ 
+ene -(R„+Rs)T c ^ _ £e|n + ee|n e~AT c ) 

1 r , p r -R„.T, n r -AT,ul- e -(^+^+ A ) T - 1 _ e -(2fl s +^+A)T c 
n l e s v n Rs+R^+X 2R S +R l _ l +X J 
i , . \E> „-i? u T c rl-e-( fls+fl M) T c i_ e -(2fl s +-R M )T Cl 
v " n l e ' A L Rs+R^i 2R S +R t _ l J 
, , d „-fA,+A)T,rl-e-C^+«M)Tc i_ e -(fls+fl M +A)T Cl 

fl, l e L Hs+H^j, Hs+ti^+A -i 
x„-R„T,rl-e-^+^+ A ) T - i_e-( 2 «s+«M+*)T Cl1 


R{s) 


Rs-start ' P(Q\RsT c ) 


R(e) 


Restart ' P(0| i2 fl T c ) ■ [1 - 6 n | e (l - e" ATc )] 


R{n) 


Rn-start ' -P(0|-R s T c ) 


R(en) 


i?e- si art ■ (1 - e" ATc ) • P(0|fl s T c ) • e n , e 


R(ss) 


Rs-start ' -P(ll-RsT'c) 


R{se) 


i? s _ start • P(0|i? s T c ) ■ e e {(l 6„| e )(l e -*« T + (e- fl ' flflT - e"^-)} 


R(es) 


i? e - s tart • P(1\R S T C ) ■ [1 - e n | e (l - e~ ATc )] 


R(ns) 


Rn-start ' P (1 -^s^c) 


R(sn) 


i? s e-^ T = • e n [1 e -(«iB D +A)T c i. 

{ R R r R (i e d j[i e -(«M+^m] + «m . ^ e - ( ij M+ fl s+ A)T c]+ 

e -(i? M +i? s )T c . ^ _ + . e -AT c ) + fi 3 . e -i? M T c . 

[1 e~ A ^ A +A -(l e -(^+2i? a+ A)T c)]} 


R(sss) 


Rs— start ' ■f > (2|-R s r c ) 


R(ess) 


i?e- s tart ■ P(2\R S T C ) ■ [1 - e n | e (l - e" ATc )] 


R(nss) 


Rn— start ' P^RgT^ 


R(ses) + R(sse) 


R(se) ■ R S T C 


R(sns) + R(ssn) 


R(sn) ■ R S T C 


R(sen) 


R s - Start ■ P(0\R s T c ) ■ e e e nle [l - e ~ R ^ - 1 R ^{ e - R '^ - e" A ^)] 


R(esn) + R(ens) 


R(en) ■ R S T C 



TABLE III: Analytic event rate expressions for all 1-fold, 2-fold and 3-fold coincident events. 
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